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^\ 1 We briefly review recent developments in black hole accretion disk theory, placing new 

emphasis on the vital role played by magnetohydrodynamic (MHD) stresses in trans- 
porting angular momentum. The apparent universality of accretion-related outflow phe- 
nomena is a strong indicator that vertical transport of angular momentum by large-scale 
MHD torques is important and may even dominate radial transport by small-scale MHD 
turbulence. This leads to an enhanced overall rate of angular momentum transport and 
allows accretion of matter to proceed at an interesting rate. Furthermore, we argue that 
when vertical transport is important, the radial structure of the accretion disk is modified 
and this affects the disk emission spectrum. We present a simple model demonstrating 
' that energetic, magnetically-driven outflows give rise to a disk spectrum that is dimmer 

' and redder than a standard accretion disk accreting at the same rate. We briefly dis- 

cuss the implications of this key result for accreting black holes in different astrophysical 
systems. 
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PACS Nos.: include PACS Nos. 



1. Introduction 

Accretion dislcs are ubiquitous in astropliysics. They form around newly born stars, 
in stellar binary systems and around supermassive black holes in the nuclei of galax- 
ies. Their eiBciency at extracting gravitational energy is truly remarkable: in the 
case of a rotating black hole, the conversion efficiency can exceed a staggering 40% 
of the equivalent rest mass energy of the accreting matter. It is now widely accepted 
that the most luminous sources in the Universe, namely quasars and gamma-ray 
bursts, must necessarily be powered by accretion onto a black hole. It is also widely 
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accepted that high-energy astrophysical sources, such as X-ray binaries and Active 
Galactic Nuclei (AGN), are also accretion-powered. 

The standard theory of astrophysical disk accretion was formulated over thirty 
years ago. Since then, arguably the most important theoretical milestone has 
been the demonstration via numerical simulations that the removal of angular mo- 
mentum required for accretion to proceed is mediated by magnetohydrodynamic 
(MHD) turbulence driven by the weak-field magnetorotational instability (MRI - 
see Ref.[4]for a review). Notwithstanding these groundbreaking developments, it still 
remains unclear precisely how energy can be channelled from the accretion flow to 
an external diffuse region where high-energy processes can operate. This outstand- 
ing issue is inextricably linked with the ubiquity of outflow phenomena associated 
with accretion disks. Indeed, the formation and nature of relativistic jets remains 
one of the most formidable theoretical problems in this field. 

In this brief review, we highlight some recent theoretical progress made towards 
understanding vertical transport in accretion disks. ^A macroscopic, mean field ap- 
proach is adopted for the MHD turbulence. This allows us to identify the dominant 
mechanisms responsible for angular momentum transport as well as the main con- 
tributions to the global energy budget of the system. We also demonstrate quanti- 
tatively how outflows, both Foynting-flux-dominated and mass-flux-dominated, can 
modify the spectrum of disk emission substantially from that predicted by standard 
accretion disk theory. 

2. MHD Disk Accretion 

In the following, we present statistically-averaged equations in a cylindrical (r, </), z) 
coordinate system for a fluid which is time-independent and axisymmetric in the 
mean {{d/dt) — {d/d(f)) — 0, where the angle brackets denote ensemble-averaged 
quantities). We use Newtonian physics throughout, with the gravitational potential 
i^G = —GM{r^ + z^)~^/^. The fluid description is valid down to the innermost 
marginally stable orbit rj « 6rg, where Tg — GM/c^ is the gravitational radius 
of a black hole of mass M. The mean-field conservation equations are integrated 
vertically over an arbitrary disk scaleheight, h = h{r). Quantities calculated at 
the disk surface (z = ±.K) are denoted by a ± superscript, X^, and we assume 
reflection symmetry about the disk midplane, so that |^^| = Midplane values 

of variables are denoted by Xq. We assume h is much less than the radius (the "thin 
disk" approximation), so that quantities of order h/r and dh/Ar are neglected. 

We adopt a mass-weighted statistical averaging approach in which all variables 
are decomposed into mean and fluctuating parts, with intensive variables such as 
the velocity mass averaged according to 

v^=v, + v[ , {pv[)^Q , (1) 

while extensive variables, such as density, pressure and magnetic fleld, averaged the 
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following way: 

P^p + p' , {p')=0 , (2) 
p = p + p' , {p')=0 , (3) 
B, = B, + B[ , {B[) = . (4) 

Note that intensive averages are denoted by a tilde, extensive averages are denoted 
by a bar and fluctuating components are denoted by a prime. The fluctuating ve- 
locity components are restricted to subsonic speeds because the MRI is a weak-field 
instability which drives subsonic turbulence, with {pv'^) < (pc^) <C pr^^l^^, where 
Cs — {kT / p.mpY/'^ is the local sound speed and JIk = (GM/r'^)^/^ is the keplerian 
angular velocity. The only restriction we place on the mean fluid velocity compo- 
nents is that they satisfy "D^ ^ Vr^Vz and that and are independent of z, 
simplifying the vertical integration. The combined magnetic and Reynolds stresses 
are defined by 

{UJ)^{^)~5^,{^)~{pvW,) (5) 

In what follows, we summarize the salient equations for radial and vertical trans- 
port of mass, angular momentum and energy. A full derivation of these equations 
can be found in our earlier paper. ^ For clarity, we have simplified the presentation 
of the relevant equations by omitting negligible correlation terms and by remov- 
ing the notation for averaged extensive and intensive quantities. Thus, averaged 
quantities are implicitly assumed. 



2.1. Mass transfer 

Vertical integration of the mean-field continuity equation gives 



— / 2TrrpVr dz + Airrp'^v^ = . (6) 



We now introduce the usual definitions for the surface mass density, 

S(r)^ r%dz (7) 

J-h 

and mass accretion rate, 

Ma(r) = 27rrE(-w,) . (8) 
We also introduce an analogous mass outflow rate, 

/•oo 

Mw(r) = / Anrp+v^dr (9) 



associated with a mean vertical velocity at the disk surface, i.e. at the base of 
a disk wind. This is the total rate at which mass leaves the disk surface at radii 
exterior to r. 
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In terms of the above definitions, the vertically integrated continuity equation, 
can be written as 

^Ma(r) = A7:rp+vt = -^^^(r) (10) 

implying that 

Ma(r) + M„(r) = M^in) + M„(ri) = constant = M , (11) 

where Ma(ri) + M„(ri) is the total mass flux at the innermost stable orbit, n. 
Equation (fTO| implies that under steady-state conditions, the radial mass inflow 
decreases towards small r at the same rate as the vertical mass outflow increases in 
order to maintain a constant nett mass flux, M , which is the nett accretion rate at 
r = oo, i.e. M = Ma(oo). 

2.2. Angular momentum 

The azimuthal component of the momentum equation is: 

1|: (r^pvrv,) + §- ipv.v.) = 1|- (r^iU,)) + ^ . (12) 



Integrating this equation over z and applying mass continuity (|10p gives 

+ 2nr^Tr 



d 
dr 



r'n^~A7rr^{U^)+ , (13) 



where 



Tr^ = I {trd,) dz (14) 
-h 



is the vertically integrated r(j) stress. The terms on the the left hand side of eqn. ([T3|) 
describe radial transport of angular momentum associated with radial inflow (ac- 
cretion) and MHD stresses acting over the disk height. The terms on the right hand 
side describe vertical transport of angular momentum resulting from mass loss in a 
wind and MHD stresses on the disk surface. The magnetic part of the MHD stresses 
are given by {U^) {BrB^)/ATT and (t^^) {B^B^)/4:Tt (c.f. eqn. E]) . 
Radially integrating equation (fT3|) (taking w J7k) gives 



dr. 



(15) 

This is a generalized conservation equation for angular momentum in accretion 
disks. It can be equivalently expressed in terms of the flux of angular momentum, 
J, as follows: 

ia(r) - Mn) = jr{r) - Mn) + Ur) , (16) 

where Ja = M^r^^l is the angular momentum flux of the accreting matter, Jr = 
—2i:r'^Tr^ is the radial flux of angular momentum carried by the rcj) stresses and = 
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r^Q dAfa/dr — 47rr^ {t(t>z)^ dr is the vertical flux of angular momentum carried 
by outflowing matter and by the (jjz stresses. Thus, both radial and vertical transport 
can contribute to the nett rate at which matter loses angular momentum and moves 
radially inwards. In other words, the mass accretion rate, Ma, is determined by the 
nett rate at which angular momentum is transported radially outwards by MHD 
stresses acting over the disk height and vertically outwards by MHD stresses acting 
over the disk surface. 

The most powerful astrophysical sources are inferred to be accreting matter 
at very high rates. Quasars, for instance, radiate energy at an astonishing rate, 
typically 10^^"'** erg s~^. Much of this luminosity can be identifled with thermal 
emission from an accretion disk around a supermassive black hole. If all the ex- 
tracted gravitational binding energy of the accreting matter is converted into disk 
radiation, then the radiative efficiency is just equal to the accretion efficiency, rj, 
which is typically 10%. Thus, the rate at which matter needs to be accreting is 
Ma = L/{ric^) ~ (0.2 — 2OO)A/0yr~^. Identifying the dominant mechanisms re- 
sponsible for such high mass accretion rates has been a longstanding challenge 
confronting accretion disk theory. This problem is made even more severe if one 
considers in addition conversion of energy into non-radiative forms, such as rela- 
tivistic jets, which can carry a total kinetic power at least as comparable to the 
radiative disk luminosity. 

It is interesting to compare the rate of angular momentum transport indicated 
by (fT5|) with that predicted by standard accretion disk theory. 1^1 In a standard 
disk, angular momentum is transported radially outwards only and the stresses 
responsible for this are parameterized by a dimensionless parameter a such that 

Tr4, = ac,/iSr— (17) 
or 

with < a < 1 for stresses due to subsonic turbulence. In addition, the stresses are 
assumed to vanish at rj and Ma is necessarily constant with r (since there is no 
vertical mass flux). Thus, the angular momentum equation, (jl5[) . reduces to 



1- ( - 

n 



^ -1/2' 



27r r acs hY, 



dr 



(18) 



There are two important limitations to the a-disk formalism. Firstly, a is treated as 
a global parameter, so that it remains constant at all radii; there is no a priori reason 
why the radial distribution of turbulent stresses should be parametrized in this way. 
Secondly, because the a— disk model remains a phenomenological prescription, it is 
unclear whether the values of a needed to account for the inferred accretion rates 
in different sources correspond to a physically plausible realization of the actual 
turbulent stresses in real accretion disks. 

Numerical simulations may be uniquely capable of addressing these issues. Sim- 
ulations of accretion disks have demonstrated conclusively that turbulent MHD 
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stresses can indeed remove angular momentum from matter, thus facilitating the ac- 
cretion process (see Ref. |4]for a review). However, 3D simulations of MRI-generated 
MHD turbul ence in accretion disks have so far been unable to produce high mass ac- 
cretion rates. 1^^^^ On the other hand, high accretion rates are recovered from 
3D MHD simul ations when a large-scale, poloidal mean magnetic field is explicitly 
included. I ^^ l ^^ l These simulations and others'^, as well as semi-analytic models 
show that at small radii, vertical transport of angular momentum by a large-scale 
magnetic torque is more efficient than radial transport by MHD turbulence. 

The numerical results indicate that MHD turbulence, whilst important for the 
microphysics of accretion disks, cannot be solely responsible for the efficacy of an- 
gular momentum transport. Indeed, as we show explicitly in eqn. IjlSp . angular 
momentum can be removed vertically from an accretion disk by the mean fluid 
and magnetic fields (i.e. nonzero Vz and {B^Bz) components at the disk surface - 
see also Ref. fT9| . We can make a direct comparison between the value of Ma pre- 
dicted by the standard model (which considers only outward radial transport) and 
that predicted by our model (which includes outward vertical transport). Using the 
a-disk formalism, the solution for Afa generalizes to 



r 

Ma = 27rQ;Cs/lS — 



dr 




1- ( - 



. 1/2' 

r ' 



(19) 



which is larger than that predicted by the standard disk model, given by (|18p . by 
a factor 1 -I- J^j Jr- That is, our model predicts that in regions where the vertical 
angular momentum flux exceeds the radial angular momentum flux, the local rate 
of mass accretion is larger than that predicted by standard accretion disk theory. 

We conclude here that it is large-scale, organized MHD stresses in the form 
of magnetic torques and mass outflows, rather than small-scale, turbulent MHD 
stresses, that lead to an enhanced overall rate of angular momentum transport (and 
hence, high Ma) in astrophysical accretion disks. The fundamental relationship be- 
tween angular momentum and mechanical energy {dE = f2dJ) then also implies 
that large-scale MHD effects must also be responsible for driving high-energy phe- 
nomena associated with accreting sources which must necessary originate outside 
the disk in a relativistic jet or magnetized atmosphere (corona). In the following 
section, we examine the global energy budget of accretion disks in which angular 
momentum transport is prescribed by the generalized conservation equation (I15|) . 



2.3. Radiative disk flux 

We now consider energy conservation in MHD disk accretion. Again, we refer the 
reader to our earlier paper (Ref. \5} for details of the derivations. 

Accretion power is the rate at which gravitational binding energy is extracted 
from the accreting matter. This energy can be converted into mechanical (e.g. ki- 
netic, Poynting flux) and non-mechanical (e.g. radiative) forms. The rate at which 
this occurs is determined by keplerian shear in the bulk flow, s^^ — ^rdfl/dr, with 
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dn/dr 



jfj/r. The rate per unit disk surface area at which energy is emitted in 



the form of electromagnetic radiation is determined by the internal energy equation. 
If there are negligible changes in the internal energy and enthalpy of the gas, then 
the disk radiative flux, Fd, is approximately equal to the rate of stochastic viscous 
dissipation of the turbulent energy, which occurs on the smallest scales, either at the 
end of a turbulent cascade or via magnetic reconnection. If there is negligible trans- 
port of turbulent energy from the source region, then turbulent energy is locally 
dissipated at a rate equivalent to its production rate, {tij)sij « {trtt>)srit> (see Ref. [5] 
for a discussion on the relative importance of production, transport and dissipation 
of turbulent energy). The internal energy equation then implies 



(20) 



Vertically integrating over the disk height yields the following expression for the 
radiative flux emerging from the disk surface: 



1 dn 



(21) 



Note that in the a-disk formalism, defined by (fT7|) . F'^ is directly proportional to a. 
In a standard disk, a is also directly proportional to Ma (c.f.[T8]). That is, a standard 
a-disk model predicts that the disk surface brightness varies linearly with the mass 
accretion rate. Our model, in contrast, predicts that disks can be accreting at a 
relatively high rate, but locally emit relatively little radiation if vertical mean-field 
transport is locally more efficient than radial turbulent transport. 

According to our model, the level of the turbulent MHD stresses T^^ available for 
internal dissipation is determined by the relative efficiency with which the extracted 
accretion energy is converted into non-radiative (mechanical) form. In other words, 
is specified by the angular momentum conservation relation (fT5|) . which gives 



■ T;0(r) 



MSI 

271 

1 

27rr2 



1 



Ma(ri) /n\i/2 



A4(r) vr 



7j r,^(n) 



dr 



dr 



(22) 



Substituting this expression into the internal energy equation (|2ip yields the fol- 
lowing solution for the disk radiative flux: 



^^d+(^) 



3GMMa(r) 



STrr-^ 



1 



Ma(ri) /ri^i/2 
Ma(r) 



(7) 



dr (23) 



87rr2 



dMa 
dr 



47rr2(t^,)- 



The first two terms on the right hand side of this equation describe the rate at which 
gravitational binding energy is extracted from matter as it accretes (i.e. loses angular 
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momentum); the second term in particular describes the rate at which nonzero MHD 
stresses at the innermost stable orbit locally dissipate turbulent energy (in practice, 
a convenient inner boundary condition is to set this term equal to zero, although 
in principle, energy can still be extracted beyond this boundary) . The last term on 
the right hand side of ([23|) describes the rate at which energy is removed through 
the disk surface by outflowing mass and Poynting fluxes. 

The result ([23|) for the radiative flux of an accretion disk corrected for the effects 
of outflows can be expressed as 



where 

Uir) = 

is the accretion factor and 
Mr) 



[Mr) -Mr)] 



A4(n) 

Ma(r) 



Ma(r)r2f7 



1/2 



.dMa 
dr 



27:rfTr4,in) 
M^{r)r^VL 



47rr2(t^^}- 



dr 



(24) 



(25) 



(26) 



is the outflow correction factor (the 'w' subscript denotes wind). Note that /w is 
equivalent to the fractional rate of vertical angular momentum transport in the 
disk. From eqns. (fT5)) and ([16]), we have 



/w(r) 



Ur) 
Mr) 



Similarly, 



and hence, 



Mr) 



Mr) 



Mr) -Mr) 



Mr) 

Mr) 



(27) 



(28) 



(29) 



Therefore, the efhcacy of disk radiant emission is entirely determined by the rel- 
ative importance of radial transport of angular momentum in the disk. If vertical 
transport of angular momentum dominates, then energy is efficiently removed from 
the disk before being locally dissipated therein. 

It is noteworthy to compare our outflow-modified disk flux with the disk flux 
predicted by standard accretion disk theory. In the absence of outflows (that is, 
when dMa/dr = and {tM~^ = 0), then /w = 0, Ma(ri) = MM) = M and hence, 
/a = 1 — (ri/r)-'^/^. In the limit of vanishing stresses at n, this is equivalent to the 
small— r correction factor from standard accretion disk theory.'' The radiative flux 
for a standard disk is then F+(r) « 3GMM^/{8TrM[l-{n/r)^/'^]. That is, standard 
accretion disk theory predicts that all the extracted gravitational binding energy is 
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locally dissipated and then radiated away. In contrast, our model predicts that the 
local internal dissipation of turbulent stresses is determined by the relative local 
rates of extraction and removal of gravitational binding energy. In regions where 
vertical transport of angular momentum dominates radial transport, the extracted 
energy flux is removed from the disk more efficiently than it is locally dissipated. 
Thus, a fraction of the nett accretion power is chanelled into MHD outflows. 

3. An Outflow-Modified Multi-Colour Disk Model 

An expression for the total disk radiative luminosity, L^, is obtained by integrating 
each term in eqn. (|23p over all disk radii, from r = n to r = oo. Similarly, the 
spectrum of disk emission can be calculated assuming local blackbody emission and 
summing the spectrum from each annulus. This is the Multi-Colour Disk (MCD) 
model used for standard disks. ^201^^3 t^[1\ henceforth refer to our generalized disk 
prescription as the Outflow-Modified Multi-Colour Disk (OMMCD) model. 

Before we can calculate the total luminosity and emission spectrum for the 
OMMCD model, it is necessary to define a specific functional form for the outflows 
from the disk surface. A simple phenomenological model for the mass accretion rate 
is (see e.g. Refs. [211 l22l [23l for other self-similar models) 

Ma(.) = ^K-)' '^-'^^ (30) 
[M ,r>r^ 

where represents a critical wind radius beyond which vertical mass loss from the 
disk surface is negligible. The wind mass loss rate then satisfies Mw{r) — M — Ma.{r), 
from the continuity equation (jlip . We can now use this model to explicitly calculate 
the factor /a, which appears as a source term in the disk flux solution (|24p and which 
is defined by ((25| . For the simplest model, we assume that the stresses Tr0(ri) at 
the inner boundary radius vanish, so that 

{/ \-(l/2+p) 
' -1/2. .-P (31) 
(^) '^^^w 

We can define a similar model for fwir), given by ([26|) . From the vertically- 
integrated, differential form of the angular momentum equation, (|13p . we have 

r^n^ - 47rr^t^^)+ = A (^Aiyn) + ^i27rr^Tr^) , (32) 

We wish to define a model in which the relative importance of vertical angular 
momentum transport decreases with increasing radius. We choose a simple power- 
law decline, such that 



- 47rr^t^,) + 



^(27rr^rr0) 



(33) 
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with w > 0. Thus, in the hmit w — > 0, radial transport of angular momentum by 
the internal MHD stresses Tr^ is negligible and the dominant mode of transport is 
via outflows of matter and Poynting flux. Inserting the relation ([33|) into (j32p and 
using (|30|) to simplify, we have for r < r„ 

-{1/2-p+w) 



dMa 

dr 



r\ \ r 



(34) 



and for r > r^, the corresponding solution is obtained by setting p — and 
dMa/dr = 0. This implies the following relation for the disk wind factor /w de- 
fined by (f26ll: 

1/2+p-tu J \nj ^ \rij 

-{l/2+p-w)-\ 



1/2+P 
1/2+p-w 

1/2 



1/2- 



-(1/2-^0' 



(35) 



r > r„ 



Eqns. (|3ip and ((35|) for /a and /w, respectively, now specify the solution for the 
OMMCD radiative flux as a function of radius, given by eqn. (1^^ . 

The total disk luminosity can now be calculated by integrating the flux over all 
disk annuli, as follows: 



La 



' , , „ , 3 GMM 
rAir) znr dr — 



M 



\ur)-U{r)] ^ (36) 



We can write the solution as 



La — Pp, — Pvj 



where 



Pa 



1 GMM ( 1 

2 rs 



2p 



1-p 



3p 



(1 + 2p) V r- 



-{i-py 



(37) 



(38) 



is the total accretion power and 
1 GMM ( l + 2p 



p 



l-p 



p{i + 2w) 

(i + 2p){\ + w) \ n 



— (l-p+io)' 



(39) 

is the total power removed by am MHD disk wind. Note that in the limit p — > 0, 
corresponding to neghgible mass loss (i.e. Ma — > constant), the accretion power 
approaches the solution from standard accretion disk theory. Pa = GAfMa/2ri, 
and the wind power is Pw Pa/(1 + w)- Thus, a substantial fraction of the total 
accretion power can be removed from the disk by a magnetic torque alone. This 
solution represents a Poynting flux dominated outflow and can be identified with 
the formation of relativistic jets that carry away very little matter, but transport a 
large amount of mechanical energy very efficiently across vast distances. 
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The emission spectrum predicted by the OMMCD model can be calculated in the 
same way as that for an MCD model (i.e. a standard accretion disk). Assuming each 
annulus in the disk radiates locally like a blackbody, Bi,, the total disk spectrum is 
calculated by summing the individual spectra from each annulus: 



La 



TrB^[T{r)] 27rr dr 



(40) 



where T{r) = [Fd{r)/a]^^'^ is the effective disk temperature of each annulus and a 
is the Stefan-Boltzmann constant. 



10" 



10^ 



10^ 



SS disk 



B torque 



z: ' B torque + mass outflow ' 



10 ' 

1 
I 

10 V 



10 

r/r. 



100 



Fig. 1. Predicted radiative energy flux emitted by an accretion disk as a function of radius r, in 
units of 3GMMa(r)/(87rr''), where Malr) is the local mass accretion rate (c.f. egn. I24| l. The solid 
line show the total radiative flux proflle for a standard Shakura-Sunyaev (SS) disk. The dotted line 
indicates the radial flux profile for a disk modified by a magnetized outfiow (with w = 0.1) with 
negligible mass loss (p = 0.1). The dashed line shows the radial flux profile for a disk modified by 
both a magnetic torque (with to = 0.1) and a mass-loaded outfiow (with p = 0.5). Mass loss is 
confined to radii within SOr; (see Sec. [3] for definitions of the model parameters). 



Figure [T] shows plots of the disk flux (r) for three different models of disk 
accretion around a black hole: 1. the standard "Shakura-Sunyaev"l2l (MCD) model; 
2. the OMMCD model with negligible mass outflow {p = 0.1) and a magnetic torque 
carrying away « 90 % of the accretion power {w = 0.1); and 3. an OMMCD model 
with a stronger mass outflow (with p — 0.5) as well as a magnetic torque (with 
w = 0.1). The mechanical energy carried away by this outflow is also « 90 % of the 
total accretion power. In all cases, mass loss is restricted to radii r < 50ri. Model 
2 corresponds to a disk modified by a Poynting flux dominated outflow such as a 
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10 




10" 10'" 10'= lO" 10" 10'' 10" 
v(Hz) 

Fig. 2. Top: Predicted emission spectra for an accretion disk around a Af = 10 Mq black hole 
accreting at a rate of M = 2.7 X 10~^ Mq yr~^ at r = oo. Bottom: Predicted spectra for M = 
IO^'Mq and M = 0.27 Mq yr-^ at r = oo. In both cases, solid lines show the spectra for a standard 
Shakura-Sunyaev (SS) disk. Dotted lines indicate the spectra for a disk modified by a magnetic 
torque, with w = 0.1 and with negligible mass loss (p = 0.1). Dashed lines indicate the spectra 
for a disk modified by both a magnetic torque (with ui = 0.1) and a mass- loaded outfiow, with 
p = 0.5. Mass loss is restricted to radii within SOn (see Sec. [3] for definitions of these parameters 
and Fig. [T] for the corresponding local radiative fluxes). 

magnetized corona and/or jet. In model 3, Pa is lower (by about a factor of 2) than 
that in the first two models because there is less matter reaching ri. 

The corresponding disk emission spectra L^^u are shown in Figure[2]for a stellar- 
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mass black hole with M = IOMq and accretion rate M — 2.7 x lO^^Afoyr^^ at 
r = oo and for a supermassive black hole with M = IO'^Mq and M — 0.27 Mq yr~^. 
These accretion rates represent the nett mass influxes at r = cxd and correspond 
to the Eddington rates when the radiative efficiency is rj — 0.1. This is strictly 
only the case for the SS disk model. When outflows are important, either due to 
a magnetic torque alone or in addition to a mass-loaded disk wind, the radiative 
efficiency can fall below the nominal 10% predicted by standard theory because the 
outflows remove energy from the disk (as well as angular momentum), so there is 
relatively less energy to dissipate and radiate away. 

As is clearly evident in Figs. [Hand [21 outflows can substantially modify the local 
disk emission profile and hence, the overall spectrum of emission. The generalized 
OMMCD model predicts a high-energy cut-off in the emission spectrum at lower 
energies than that predicted by the standard SS model. This results from the im- 
portance of outflows at small radii, where the highest-energy emission originates. In 
the case of stellar-mass black holes, the spectrum is most affected at X-ray energies 
(between 10"'^'^"^^ Hz), while for supermassive black holes, it is the extreme ultravi- 
olet region (^ 10^^ Hz) that is most affected by outflows. The OMMCD model also 
predicts a broadband region of the disk spectrum (around visual frequencies) that 
is much flatter than the characteristic i'^^^ law predicted by standard theory. 

Our results demonstrate that the modifications to a standard disk spectrum as 
a result of energetic outflow phenomena are clearly not negligible when the outflows 
are primarily responsible for the removal of angular momentum from accreting 
matter. Direct observational verification of the predicted OMMCD spectrum may 
be possible for only some types of accreting sources, however, as the outfiows will 
also produce their own characteristic radiative signatures that may overlap in the 
spectral energy band where we predict the disk to be most strongly modified. This 
can confuse our interpretation of the emergent observed spectrum. 

In the case of galactic X-ray binaries (XRBs), for instance, the disk spectrum can 
be substantially modified at X-ray energies (typically above 1 keV, corresponding 
to frequencies above 10^^ Hz - see Fig. [2]). These sources often exhibit a power-law 
X-ray spectrum above 1 keV that is thought to arise in a magnetized, tenuous atmo- 
sphere (a "corona") and/or relativistic jet as a result of inverse Compton scattering 
off disk photons (see Ref . [24l for a review). Existing models for this emission process 
simply use a standard disk spectrum for the seed photon distribution. However, this 
is clearly not a self-consistent calculation since the formation of a corona and/or jet 
from disk magnetic fields will inevitably modify the disk radial structure and hence, 
photon spectrum, as we have shown here. Thus, although we would not expect to 
directly see the modified disk spectrum in XRBs, our model can be indirectly tested 
by using the underlying OMMCD spectrum as the source of seed photons that are 
upscattered in a disk corona and/or jet. By comparing the predicted and observed 
X-ray spectra, we can improve our interpretation of the source characteristics and 
obtain tighter constraints on key physical parameters such as mass accretion rate. 
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The above type of c alculation has been performed by us for ultra-luminous X- 
ray sources (ULXs). These are exceptionally luminous XRBs found in external 
galaxies. They probably involve accretion onto a black hole more massive than 
those found in galactic XRBs. They provide an effective test for our modified disk 
model because some ULXs exhibit a low energy ( < 1 keV) spectral component that 
can be interpreted as disk emission. However, because these sources emit virtually 
all of their radiative power as hard, power-law X-rays, which must necessarily be 
produced outside the disk, the accretion disk must be substantially modified and a 
standard disk model should not be used to fit the soft spectral component. Indeed, 
we show - that the OMMCD spectrum can adequately fit the observed soft spectral 
component, implying a black hole mass M ~ IOOMq. 

A similar test can also be performed for quasars and other AGN, which are 
powered by accretion onto a supermassive (10^^^ M0) black hole. These sources 
exhibit a prominent optical-ultraviolet (UV) continuum feature known as the "big 
blue bump". ^' This e nerget ically significant feature is generally interpreted as 
accretion disk emission. I ^^ l ^^ l However, attempts at fitting accretion disk spectral 
models to the observed spectra have had mixed success (see Ref.ISOlfor a review). For 
AGN, the disk spectrum will be most strongly modified by outflows at ultraviolet 
(UV) wavelengths (see Fig. [5]). Although there is less overlap between the modified 
disk and corona/jet X-ray emission than in the XRB case, the poor transmission 
of UV radiation through Earth's atmosphere means that direct observations of 
the predicted OMMCD spectrum in AGN are not straightforward. Nevertheless, 
spectral energy distributions have been compiled from optical observations of high- 
redshift quasars and from satellite UV observations with the Hubble Space Telescope 
(HST) and the Far Ultraviolet Spectroscopic Explorer (FUSE). Interestingly, these 
observations have revealed a far-UV break in the spectral energy distributions of 
quasars (see Ref. [31] and references cited therein) . Specifically, the observed spectra 
decline dramatically at wavelengths shorter than 1100 A (i.e. frequencies above 
3 X 10^^ Hz). This is difficult to explain with standard accretion disk models, which 
predict that the spectrum should continue to rise into the extreme UV (see Fig. [2]) , 
but is consistent with the predictions our modified disk model (see Ref. l32| . 

4. Concluding Remarks 

The standard theory for astrophysical disk accretion has enjoyed remarkable suc- 
cess since its inception over 30 years ago. However, ongoing rapid advances in in- 
strumental technology are placing increasingly tighter observational constraints on 
theoretical models that use a standard accretion disk. The ubiquity of energetic 
outflow phenomena in accreting astrophysical sources gives us two important clues 
to the underlying processes responsible for accretion: 1. vertical transport of angular 
momentum cannot be ignored and may indeed be the dominant mode of transport 
at small radii; and 2. the removal of accretion energy from the disk must result in a 
modified radial disk structure, particularly at small radii, and consequently, a disk 
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spectrum that is dimmer and redder than a standard disk accreting at the same 
rate. We have presented a simple model that can be used to quantitatively calculate 

model spectra for disks modified by a magnetized jet and/or corona and by mass- 
loaded winds. These model spectra can be applied to a variety of different sources 
to improve our interpetation of the observations and more accurately determine key 
physical parameters such as mass accretion rate and black hole mass, as well as the 
partitioning of accretion power into radiative and non-radiative forms. 
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